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Abstract. We have made self-consistent models of the density and temperature profiles of the gas and dust surrounding em- 
bedded luminous objects using a detailed radiative transfer model together with observations of the spectral energy distribution 
of hot molecular cores. Using these profiles we have investigated the hot core chemistry which results when grain mantles are 
evaporated, taking into account the different binding energies of the mantle molecules, as well a model in which we assume 
that all molecules are embedded in water ice and have a common binding energy. We find that most of the resulting column 
densities are consistent with those observed toward the hot core G34.3+0.15 at a time around 10"* years after central luminous 
star formation. We have also investigated the dependence of the chemical structure on the density profile which suggests an 
observational possibility of constraining density profiles from determination of the source sizes of line emission from desorbed 
molecules. 
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1. Introduction 

Hot cores are compact (diameters < 0.1 pc), dense (n> 10^ 
cm"^), hot (T> 100 K), and dark (A^ > 100 mag) molecular 
cloud cores. At least some of them are believed to be precursors 
to ultracompact HII regions and in the early stage of forming 
massive stars because luminous infrared sources, massive bipo- 
lar outflows and/or inflow motion are observed associated with 
them (e.g., Henning et al. 2000; Kurtz et al. 2000; Churchwell 
2002). They are observationally characterized by anomalously 
large abundances of some molecular species including hydro- 
genated, saturated, small molecules such as H2O and H2S, and 
large complex molecules such as HCOOCH3 and (CH3)20, in 
contrast to colder dark clouds. It is thought that these molecules 
are produced under conditions specific to star formation: the 
freeze-out of gas-phase species onto dust grains and the grain 
surface reactions in the cold collapsing pre-protostellar phase, 
and the evaporation of the grain mantle material and the sub- 
sequent gas-phase reactions in the hot protostellar phase (e.g., 
Millar 1993, 1997). Therefore, the observed large column den- 
sities of some molecules in hot cores are also expected to be an 
indicator of luminous high mass star formation which yields a 
large mass of high temperature gas. Previously, chemical mod- 
els of hot cores have been studied together with the pre-stellar 
phase or with some initial conditions of evaporated grain man- 
tle species, and have been successful in reproducing their ob- 
servational properties (e.g.. Brown et al. 1988; Millar et al. 
1991; Charnley et al. 1992; Caselli et al. 1993; Charnley 1997; 
Miflar et al. 1997; Viti & Wifliams 1999; Rodgers & Charnley 
2001; Doty et al. 2002). 



Physical properties of hot cores have been investigated on 
the basis of observational studies of molecular line and dust 
continuum emission. Some physical models together with ra- 
diative transfer calculation and/or chemical calculation have 
succeeded in reproducing the observations of spectral energy 
distributions and/or molecular column densities (e.g., Millar et 
al. 1997; Kaufman etal. 1998;Osorioetal. 1999; Hatchell et al. 
2000; van der Tak et al. 2000a; Doty et al. 2002). On the other 
hand, the dynamical processes and the formation mechanisms 
of massive stars are not yet established because of the diffi- 
culty of treating in a self-consistent manner young stellar evo- 
lution, the accretion flow and radiation pressure. Non-spherical 
accretion to protostellar cores, for example, through accretion 
disks, and/or coalescence of lower mass stars may be involved 
in massive star formation, but there are, as yet, no accepted 
theories (e.g., Bonnell et al. 1998; Stabler et al. 2000; Yorke & 
Sonnhalter 2002). The nature and timescale of the dissipation 
of the surrounding envelopes of young massive stars, for ex- 
ample, by powerful outflows and ionizing strong winds which 
stop the accretion, are also uncertain but important in deter- 
mining the mass of forming stars in the case of the isolated star 
formation (e.g., Shu et al. 1987; Nakano et al. 1995). 

In this paper, we have modeled the density and tempera- 
ture profiles of hot cores, making use of radiative transfer cal- 
culations, and then investigated the hot core chemistry, taking 
into account the temperature dependent grain mantle evapora- 
tion. In the following section, we present the physical model of 
the density profiles constrained by observational results and the 
temperature profiles calculated by solving the radiative trans- 
fer equation. In Sect. 3 we introduce the chemical model along 
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with the temperature dependent grain mantle evaporation pro- 
cess. These models are applied to the hot core G34.3+0.15, 
and the results of chemical calculation are described and com- 
pared with the observations in Sect. 4. Also, we investigate the 
dependence of the chemical structure of hot cores on (1) the 
trapping of mantle molecules in water ice and (2) the density 
profile generally to discuss the chemical and physical condition 
of massive-star-forming clumps in Sect. 5. Finally, the results 
are summarized in Sect. 6. 

2. Physical model 

In this paper we consider a spherically syrmnetric molecular 
clump illuminated by a central star with luminosity L,. The 
density profiles modeled are constrained by observational re- 
sults (Sect. 2.1), and the temperature profiles are calculated 
self-consistently under the assumption of radiative equilibrium 
by solving the radiative transfer equation (Sect. 2.2). 

2.1. Density profiles 

As we will see in Sect. 5, the density profile of a clump is 
one of the crucial properties which afi"ects its chemical struc- 
ture. However, there is no established theoretical model be- 
cause the dynamical process of massive star formation is still 
unknown due to the difficulty in comprehensive treatments of 
young stellar evolution and accreting flow as well as ionizing 
radiation, radiation pressure, stellar wind, bipolar outflow, and 
so forth (e.g.. Stabler et al. 2000; Yorke & Sonnhalter 2002). 
Thus, in this paper we model the density profiles of clumps 
simply based on the submilhmetre observations by Hatchell et 
al. (2000), avoiding uncertainties in theoretical models. 

Hatchell et al. (2000) investigated spectral energy distribu- 
tions (SEDs) and radial profiles of 450 and 850 /xm contin- 
uum dust emission of five ultracompact HIT regions to conclude 
that chemically rich hot molecular clumps can be modeled by 
a combination of compact high density cores plus r den- 
sity profile envelopes. Also, van der Tak et al. (2000a) stud- 
ied 14 envelopes around young massive stars by means of 
CO, CS, and H2CO line emission maps, and SEDs and ra- 
dial profiles of submm dust emission to find that the best-fit 
values of the power law indices for density profiles of the en- 
velopes range from -1.0 to -1.5. Theoretically, an isothermal 
self-gravitating gas is known to have r"-^ density profile via the 
balance between self-gravity and pressure gradient forces (e.g., 
Ebert 1955; Bonnor 1956; Larson 1969; Shu 1977). If the non- 
thermal velocity dispersion, which has been observed in Une 
widths towards various sizes of interstellar clouds and known 
to have a scahng-law of Sv oc r"^ (e.g., Larson 1981; Falgarone 
et al. 1992; Caselli & Myers 1995), affects the pressure gra- 
dient force, it may flatten the radial density profiles of clumps 
hke the logotrope model (e.g., McLaughhn & Pudritz 1996). 

Taking into account of these observational and theoretical 
studies, we use the following spherically symmetric hydrogen 
number density distribution in this paper: 



where the radius of ro- = 0.05 pc is adopted for the transition 
between isothermal and nonthermal regions, according to the 
observations of molecular line widths (e.g.. Fuller & Myers 
1992; Myers & Fuller 1992). The hydrogen column density of 
the clumps, A^h = 2 J^°°'n(r)<ir, is fixed as A'h = 10^^ cm"^, 
based on the observations of the chemically rich hot molecu- 
lar clumps by Hatchell et al. (2000). The inner radius is set at 
the stellar surface, /?„ and the outer radius is set at rout = 1.0 
pc. The central densities «o and ni are determined to satisfy the 
relation of no/(l + ''tr/rc)^ = «i/(l + rtr/rtr)'^ and the column 
density of A^h = 10^^ cm'^. 

2.2. Temperature profiles 

The temperature profiles are calculated self-consistently, as- 
suming that the heating source of the clump is a central star 
with luminosity L„ under the condition of local radiative equi- 
librium, that is, the balance between absorption and reemission 
of radiation by dust grains at each point in the clump: 



4;r 



dvKyBy[Tir)] = I dvKy(T)dnd<pIyir,n,^), (2) 
Jo J 

where T, Ky, and By(T) represent the temperature, the 
monochromatic opacity, and the Planck function for black- 
body radiation at a frequency v, respectively. Local thermody- 
namic equilibrium, rjy = KyByiT), is adopted, where rjy is the 
monochromatic emissivity. The specific intensity ly is calcu- 
lated by solving the radiative transfer equation. 



Iy(r,i^,cp)= Ky(r')p(r')By[T(r')]e 
/() 



-Ty(r' 



(3) 



where Ty(r') is the specific optical depth from r' to r. The short 
characteristic method in spherical coordinates (DuUemond & 
TuroUa 2000) is used to compute the integration. The radia- 
tive transfer equation is solved only in the dusty region (where 
the dust grains are not thermally destroyed) in this calculation. 
The inner radius for the numerical calculation is determined 
as fin = {TinlTtY^Rt, where T, and R, are the temperature 
and radius of the central star, assuming that the radial profile 
of the temperature in the dust-free region (r < ^n) is given 
by T{r) oc r which represents the profile in the optically 
thin region with gray opacity heated by a point source. For 
the dust opacity model, we adopt the frequency dependent ab- 
sorption coeflicient Ky of Adams & Shu (1985, 1986), which 
approximates the interstellar dust grains model by Draine & 
Lee (1984). The dust destruction temperature is set as T^a - 
2300K in this model (see Nomura 2002 for the details of radia- 
tive transfer calculation). In this paper the gas temperature is 
assumed to be the same as the dust temperature, based on the 
result of our calculation that the difference between the gas and 
the dust temperatures is small, except for the very inner region 
of the clump, where the gas temperature is > 10^ K due to pho- 
toelectric heating from small grains induced by the ultraviolet 
radiation from the central star (cf. Doty & Neufeld 1997). 

3. Chemical model 



= I «o/(l + r/r^f+m/il + rlr^y-^ if < 
\ni/{l + r/rcY'^ otherwise. 



It is generally believed that the relatively large abundances of 
'■^^ hydrogenated, saturated, smaU molecular species, such as H2O 
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and H2S, and large complex molecules, such as HCOOCH3 
and (CH3)20, observed in hot cores are related to the evolu- 
tion of young stellar objects and the evaporation of grain man- 
tle material. In the cold, pre-stellar phase of the dense molec- 
ular clumps the accretion timescale of gas-phase atoms and 
molecules onto grain surface, ~ 3xlO''/n [cm^^] yrat T =10K, 
is short enough since the density, n, is high via self-gravitating 
contraction of the clumps. Thus, the accretion of gas-phase 
species onto grain surfaces and subsequent surface chemistry 
produce grain mantles rich in saturated species. Once massive 
star formation occurs in the clumps, the molecules are evapo- 
rated from the grain surface and the subsequent gas-phase reac- 
tions produce the large complex molecules seen in the hot gas 
(see e.g., reviews by Millar 1993, 1997 and references therein), 
although there is some evidence that even the most complex of 
molecules may have been formed on grains. On the other hand, 
the recent detailed analyses of observations of ice absorption 
features in infrared spectra towards massive young stellar ob- 
jects have suggested that icy grain mantles in envelopes sur- 
rounding the central stars have different composition at differ- 
ent temperatures along the line-of-sight, that is, each icy man- 
tle molecule sublimates from dust grains at a specific temper- 
ature (e.g., Ehrenfreund et al. 1998). In this paper, we set the 
molecular species evaporated from grain mantles as an initial 
condition, and then calculate the subsequent time-dependent 
gas-phase chemical reactions, taking into account the different 
binding energies of the mantle molecules. 

The initial fractional abundances of species with respect to 
total hydrogen nuclei listed in Table 1 are used in this paper. 
The initial gas is assumed to be mostly molecular and slightly 
ionized by cosmic rays. The initial mantle composition is deter- 
mined so as to be consistent with the recent observations of ice 
absorption features in infrared spectra towards massive young 
stellar objects (Gerakines et al. 1999; Keane et al. 2001) and 
other hot core models (Charnley et al. 1992; Millar et al. 1997; 
Rodgers & Charnley 2001; Doty et al. 2002), and fixed so as 
to best fit the results of model calculations to observations to- 
wards G34. 3-1-0. 15 (see Sect. 4). Also, the total abundances of 
carbon, oxygen, and nitrogen are set to agree with the mean 
interstellar values of C/H= 1.4 x 10 ^*, 0/H= 3.19 x 10 ^ and 
N/H= 7.5x10-^ (Cardellietal. 1996; Meyer et al. 1997; Meyer 
et al. 1998). 

Each surface species, is expected to evaporate thermally 
with the timescale of 

Tevap,/ = Vgj. exp{Eb,i/kT), (4) 

where k and T are the Boltzmann constant and the tempera- 
ture, respectively. The vibrational frequency vo,, is given by 
vo,i - (InsEhj/n^miy^^, where «s ^ (7A)-' is the number of 
sites per unit surface area and m, the mass of the species (e.g., 
Tielens & Allamandola 1987). The binding energies E^j are 
taken from Hasegawa & Herbst (1993), Aikawa et al. (1997), 
Willacy et al. (1998), and Eraser et al. (2001) (see also ref- 
erences therein). Now, the temperature in the clump is a de- 
creasing function of distance from the central star and as- 
sumed to be unchanged with time, while the evaporation rate 
is very sensitive to the temperature as can be seen from Eq. 
0. For example, in the case of water, whose binding energy is 
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Table 1. Initial abundances and injection radii of mantle 
molecules in the model of G34. 3-1-0. 15. 



Species Abundance 


'■injj[pc] 


Species 


Abundance 


'■mj.,[pc] 


H+ 


1.0 10"" 




CO 


1.3 10-" 




He+ 


2.5 10-'2 




CO2 


3.0 10-'' 


0.18 




1.0 10-" 




H2CO 


2.0 10-* 


0.3 


Fe+ 


2.4 lO-** 




CH3OH 


2.0 10-^ 


0.12 


He 


0.1 




C2H5OH 


5.0 10-" 


0.14 


S 


5.0 10-" 




O2 


1.0 10-* 


0.5 


Si 


3.6 10-'* 




H2O 


2.8 10-* 


0.08 


C2H2 


5.0 10-^ 


0.22 


N2 


3.7 10-5 




CH4 


2.0 lO-'' 


0.7 


NH3 


6.0 10-^ 


0.17 


C2H4 


5.0 10-" 


0.21 


H2S 


1.0 10-^ 


0.29 




5.0 10-" 


0.20 


ocs 


5.0 10-'* 


0.17 



Eh,Hio/k - 5773 K, a temperature rise from T - 95K to 105K 
shortens the evaporation time by about 1/100. Thus, in this pa- 
per the evaporation process is simply modeled:- each species, 
evaporates instantaneously inside a radius rinj,, where the evap- 
oration time Tevap./ is shorter than the accretion time, 

Tacc,; = {STTa^dgnviY^, (5) 

while it is retained on the dust otherwise. That is, the gas-phase 
abundances of mantle molecules are initially set as those listed 
in Table 1 inside the injection radii rinj,,-, and zero outside these 
radii. Beyond a radius of 1 pc, the density becomes so low that 
the accretion time becomes comparable to or greater than the 
chemical time-scales so that mantle formation is not efficient. 
We do not model this extended envelope here. In Eq. (|5} n, 
a, dg, and v, - {2kT jmi)^^'^ are the number density of hydro- 
gen nuclei, the grain radius, the ratio of the number density of 
grains to n, and the thermal velocity of the species, /, respec- 
tively. The sticking probability S is set as 5 = 0.3. Eor the grain 
surface area per hydrogen nucleon, ndgO^, we used the value of 
< dga"^ >- 2.2 X 10-^^cm^ (see Rawlings et al. 1992 for the 
parameter values). The corresponding injection radii rmj,,- in the 
model of G34. 3-1-0. 15 (see next section) are also listed in Table 
1 . The binding energies of CO and N2 are low enough and the 
kinetic temperature high enough that they exist in the gas-phase 
throughout the region modelled here. 

The subsequent gas-phase chemistry is calculated time- 
dependently by means of the chemical network based on that 
of Millar et al. (1997), which consists of 209 species connected 
by 2191 reactions. The standard interstellar cosmic-ray ioniza- 
tion rate, ^ = 1.3x lO-'^s-', is adopted. The adsorption process 
of neutral atoms and molecules / onto grain surface is included 
only if Tacci < T"evap,i in this calculation. The accretion rate is 
given by 

k-dcc,i — -S na^dgnvirii, (6) 

where n, is the number density of species i. This process affects 
the gas-phase abundances slightly in the outer region, but is not 
very significant because the accretion rate of each species is 
comparable to the chemical reaction rates. 
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4. Model of hot core G34.3+0.15 

In this section, we model the physical and chemical structure of 
the hot molecular core G34.3+0.15 using the radiative transfer 
and the pseudo-time-dependent gas-phase chemistry calcula- 
tions described in the previous sections together with observa- 
tions of the spectral energy distribution and the molecular line 
emissions toward G34. 3-1-0. 15. The hot core associated with 
this source consists of a small clump of hot gas situated at 
the nose of a cometary-shaped, ultra-compact Hii region. It has 
been well-observed, particularly through spectral line surveys 
(Macdonald et al. 1996; Kim et al. 2000) at 1, 2 and 3 mm, and 
has a rich molecular line spectrum, with large abundances of 
complex, highly- saturated molecules. 



4. 1. Density and temperature profiles 

The density and temperature profiles of the hot core 
G34.3+0.15 are modeled so that the theoretical calculation of 
the spectral energy distribution (SED) is consistent with that 
observed. We take rc defined in Sect. 2.1 as a parameter and 
obtain its best-fit value by iteratively calculating the temper- 
ature profile (see Sect. 2.2) and comparing the resulting SED 
with observation. The luminosity of the central star is set as 
L, = 4.7 X lO^Lo, referring to Watt & Mundy (1999). For the 
temperature and radius of the star we adopt = 4.5 x lO'^K 
and R, = I IRq, assuming a massive zero-age main-sequence 
star (e.g., Thompson 1984; see also references therein). The 
distance to G34. 3-1-0. 15 of 3.7 kpc and the observational data 
for the SED are taken from Chini et al. (1987). As a result of 
the calculations, we find that the model with = 0.05 pc repre- 
sents the best fit to the observational data at wavelengths from 
far-infrared to millimeter. The resulting density and tempera- 
ture profiles, and the calculated SED are shown in Figs. la, lb, 
and Ic, respectively, with solid lines. The crosses in Fig.lc are 
the observational data. 

The data at wavelengths from near- to mid-infrared can not 
be reproduced by this model probably because of the inappro- 
priate assumption of the spherically symmetric density profile. 
For example, the optical depth of non-spherical cavities due 
to bipolar outflows would be less than that of the dense en- 
velopes so that if we were to observe the inner region of the 
clumps through optically thin cavities, namely, there would be 
more flux in the infrared (cf. Hatchell et al. 2000; van der Tak 
et al. 2000a). Actually, the model with the column density of 
A'h = lO^'* cm"^ (1/10 of the above envelope model) and the 
core radius of = 0.05 pc, which is plotted with a dotted 
line in Fig. Ic, better represents the observational data in the 
mid-infrared. Here, the bolometric luminosity ratio contributed 
from the envelope region with Nn = 10^^ cm"^ to the cavity 
region with A^h = lO^"* cm"^ is assumed to be 20 : 1 in order to 
fit the observational data. This will mean that the solid angle of 
the cavity will be much less than 4n and/or the radiation flux in 
the clump around the central star is not spherically symmetric 
in fact. Its density and temperature profiles are also shown with 
dotted lines in Figs, la and lb. In the following subsections, we 



:(a) 



10000 I 

0.001 



r[pcl 




Fig. 1. (a) Density and (b) temperature profiles and (c) SED 
of the model of G34.3+0.15. Combination of a envelope with 
A^H = 10^^ cm"^ (solid lines) and a cavity with A^h 



(dotted lines), whose core radii are 
SED observational data (crosses). 



: 0.05 pc, reproduces the 



only use the model with the column density ofN= 10^^ cm 



because most of the mass of the clump will be contained in this 
dense envelope region rather than the less dense, cavity region. 
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4.2. Results of chemical calculation 

4.2.1 . Time evolution of molecular column densities 

Making use of the density and temperature structures obtained 
in the previous subsection and the chemical model in Sect. 3, 
we simulate the time evolution of the hot core G34.3+0.15. 
The resulting time evolution of molecular column densities are 
shown in Fig.2. In order to compare these with the observations 




10000 100000 

t[yrl 





1e+16 


C2H5OH 




1e+15 






1e+14 


HCOOCHj 


(t) [cm ' 


1e+13 


r^'onjCcT 




1e+12 


CH^O 




1e+11 






1e+10 






1e+09 






Fig.2. The time evolution of column densities of some 
molecules in the model of G34.3+0.15. Parent molecules are 
eventually destroyed to produce daughter species around the 
timescale of 10* yr. 

of beam-averaged coluimi densities, each atomic or molecular 
column density is calculated as 



Ni(t) = 



1 



f 

Jo 



2nrNi{r, t) exp - 



\dr, 



(7) 



where Niir, f) — n(s)xi(s, t)ds is the column density of 
species ; along a line of sight at a distance r from the clump 
centre. Here, n(s) is the number density of hydrogen nuclei, 
Xiis, t) the fractional abundance of species i at time f, and s the 
depth along the line of sight. The scale size ©beam is simply 
taken as 0.25 pc which corresponds to the angular size of 13" 
at the distance of 3.7kpc, based on the half power beam width 
of the James Clerk Maxwell Telescope (JCMT) at 345 GHz (cf . 
Thompson et al. 1999). While the data of most molecular lines 
we use in this paper are from single radio telescope observa- 
tions such as the JCMT, the NH3 data is from the Very Large 
Array observation whose synthesized beam size is about 1"3 
(Heaton et al. 1989). Thus, we calculate the NH3 column den- 
sity with ©beam = 0.025 pc as well as with ©beam = 0.25 pc for 
comparison. 

From Fig.2 we can see that parent molecules, such as H2O, 
NH3, H2S, and CH3OH, which are expected to be formed on 
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CO 
C2H 

CH2C2H 

H2O 
H2CO 
HCO'' 
HCO 
NO 
NH3 
HNC 
HCN 
ON 
HC3N 
CH3CN 





"'1 1 1 '1 1 


"1 1 1 ' 






i»- 1 - 




)-#■ — I*- 














♦ 1 




4- — 1 — 






H4- — i*- 






1 






4- — if — 1 








--HK--t*- - 










♦ 1 1 






1— ♦— i.- 






1 1 






1 (♦--i^ 

...1 1 1 1 1 


■■1 1 1 



Table 2. Molecular column densities: centre vs. off-centre. 



10''^ 10^^ 



^014 .,015 ^q16 ^q17 ^q18 ^o19 ^^20 .,o21 



SO 
SO2 
CS 
OCS 
HjCS 
HCS"" 
NS 
CH3OH 
C2H5OH 
CH2CO 
CH3CHO 
HCOOCH3 
CH3OCH3 



10 



12 



10 



13 



10 



14 



N 



10^5 _ 
mol [cm"' 



10 



16 



10 



17 



10' 



Fig. 3. Comparison between the calculated molecular column 
densities from 3 x 10^ yr to 3 x 10"* yr (solid lines) and the ob- 
servations towards G34.3+0.15 (diamonds plus dotted arrows 
if they are the lower or upper limits). The results are in good 
agreement with most of the observations. 



grain surfaces, are relatively stable for timescales of at least 
10"* yr in this model. However, they are eventually destroyed 
by ionized species or atomic hydrogen to produce daughter 
species (see Sect. 4.2.3 for the details) as has been suggested in 
the previous studies (e.g., review by Millar 1997). Thus, the ob- 
served large abundances of both parent and daughter molecules 
will mean that the ages of the chemically rich hot cores are 
about lO"* yr. In particular, we compare the results of our calcu- 
lation with observations in Fig. 3, which shows that the calcu- 
lated molecular column densities from 3 x 10^ yr to 3 x 10"* yr 
(solid hues) are in good agreement with most of those observed 
toward the hot core G34.3 -1-0.15 (diamonds plus dotted arrows 
if they are the lower or upper limits), consistent with the result 
of Millar et al. (1997). Lower hmits in this Table have been de- 
rived from optically thick transitions. The calculated molecular 
column densities, with the exception of CH3CCH, HNC, H2CS, 
and CH3OCH3, satisfy the observational lower Umits or fit the 
observations to within a factor 4 at most. Discrepancies may be 
resolved in a number of ways. Dimethyl ether, CH3OCH3, may 
in fact be a mantle molecule, while the gas-phase chemistry of 
H2CS is not at all well understood. Unlike H2CO, for which 





Molecular column density [cm ^] 




Centre 


OfF-centre 


Species 


Model" 


Observation 


ModeP 


Observation 


H2 


3.6 X lO^-* 




2.0 X 10^* 




CO 


4.9 X lO^o 


> 2.9 X 10" 


2.6 X 102" 


> 1.3 X 10'" 


C2H 


7.4 X 10'5 


> 8.2 X 10'^ 


4.9 X 10" 


> 2.9 X 10'* 


CH3C2H 


1.1x10'* 


1.8 X 10'* 


1.9 X 10'3 




H2O 


3.5 X 10" 


1.3 X 10'" 


4.5 X 10'* 




H2CO 


1.7 X 10'* 


> 2.3 X 10'* 


3.0 X 10'^ 


> 1.3 X 10'* 


HCO+ 


8.7 X lO'"* 


> 1.4 X W 


6.0 X 10'* 


> 2.3 X 10'^ 


HCO 


6.4 X 10'* 


> 5.3 X 10'3 


1.3 X 10'* 




NO 


3.9 X 10'« 


> 5.8 X 10'5 


1.2 X 10'* 




NH3 


3.7 X 10'^ 




3.9 X 10'* 




NH3'' 


2.6 X 10'** 


> 2.7 X 10'* 






HNC 


6.2 X 10'^ 


< 3.3 X 10'* 


8.6 X 10'* 




HCN 


1.0 X 10'* 


1.2x 10'5 


1.3 X 10" 


> 2.3 X 10" 


CN 


5.5 X lO'* 


> 5.9 X 10'* 


3.9 X 10'* 


> 9.5 X 10'^ 


HC3N 


1.0 X 10'* 


> 2.5 X 10" 


9.7 X 10'* 




CH3CN 


5.5 X 10'3 


> 2.4 X 10'* 


6.5 X 10'2 




H2S 


1.2 X 10''' 


> 1.1 X 10'* 


1.9 X 10'* 




SO 


1.7 X 10'* 


> 24 X 10" 


2.3 X 10" 




S02 


9.3 X 10'^ 


9.3 X 10'^ 


9.9 X 10'* 


> 2.0 X 10'* 


CS 


7.5 X 10'' 


1.0 X 10'* 


14 X 10" 


> 5.3 X 10'^ 


OCS 


3.6 X 10'* 


> 1.8 X 10'* 


3.4 X 10" 




H2CS 


2.8 X 10'* 


> 3.8 X 10'* 


7.9 X 10'3 




HCS+ 


7.3 X 10'2 


> 3.0 X 10" 


3.6 X 10'^ 




NS 


4.3 X 10'* 


>4.0x 10" 


8.7 X 10'3 




CH3OH 


7.7 X 10'* 


2.6 X 10'* 


6.9 X 10" 


> 3.6 X 10'* 


C2H5OH 


2.4 X 10'5 


1.7 X 10'5 


2.3 X 10'* 




CH2CO 


3.0 X 10'* 


6.7 X 10'* 


1.3 X 10'* 




CH3CHO 


1.4 X 10'* 


24 X 10'* 


1.5 X 10" 




HCOOCH3 


1.1 X 10'^ 


1.4 X 10" 


9.2 X 10'^ 




CH3OCH3 


5.6 X 10'3 


3.7 X 10" 


5.0 X 10'^ 





^ Calculated results at 10* yr. 
In the case of ©beam = 0.025 pc. 



efficient gas-phase synthetic pathways have been studied in 
the laboratory, the analogous reactions which might form thio- 
formaldehyde do not occur. The chemistry of HNC at low tem- 
peratures is well understood; the discrepancy between model 
and observation may imply that additional high-temperature 
loss mechanisms which are not included in the model are at 
work. Finally, the methyl acetylene, CH3CCH, column density 
is about a factor of 100 less than that observed by Macdonald et 
al. (1996) from their 330-360 GHz spectral line scan. However, 
it is much closer to the value of 8.2 x lO'"* cm"^ determined 
from the 2 and 3 mm line surveys of Kim et al. (2000), similar 
to the value derived by Hatchell et al.(1998a) 20" offset from 
the hot core. Finally, we note that the formation of C3H5 , the 
precursor ion to CH3CCH, involves ions and neutrals related to 
CH4, C2H2 and C2H4. An increase in the abundances of each 
of these parents by a factor of a few would also help resolve 
the discrepancy noted above. The observational data are taken 
from MiUar et al. (1997) (see also Macdonald et al. 1996 and 
Hatchell et al. 1998a) and Ikeda et al. (2001). 
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4.2.2. Radial profiles of molecular abundances 

In order to see the effects of the temperature dependent injec- 
tion of molecules from dust grains on hot core chemistry, we 
present the radial profiles of molecular abundances at 10'' yr in 
Fig.4. We can see from the figures that the abundances of many 
parent molecules injected from grain mantles, such as C2H2, 
CO2, H2O, H2S, and CH3OH, change dramatically at ri„j„ in- 
side which the molecule i can be easily evaporated. Inside n^jj, 
most of the parents are difficult to destroy. Nevertheless, given 
their large initial abundances, even a slow destruction rate is 
large enough to produce appreciable abundances of daugh- 



(') ^ — 


CH30H 1 

- C2H50H 


HCOOCH, 
CH.CHO 
CH,DCH, 

CH2Cd - 


HCOOCH3 




. CH3CH0 / 

















r[pcl 



Fig.4. The radial profiles of some molecular abundances of 
the model of G34.3 -1-0.15 at lO'* yr. Some parent, daughter and 
other related species have dramatic changes at radii related to 
the grain mantle evaporation. 



ter products. Thus, the radial abundance profiles of the par- 
ent molecules influence those of daughter and other related 
species (see Sect. 4.2.3 for details). Abundance changes of 
some molecules along the radial direction in the clumps around 
young massive stars have been suggested by observational 
studies. For instance, the observations of the CH3CN rotational 
transition Unes indicate that they are emitted from the inner 
compact hot region of the clumps, while the observed CH3OH 
lines suggest that the gas in some clumps consists of at least 
two components, one of which is hot and abundant in CH3OH, 
and the other is cooler and less abundant (Hatchell et al. 1998a; 
MiUar & Hatchell 1998; van der Tak et al. 2000b). In addition, 
the observational abundances of the vibration-rotation lines of 
C2H2, HCN, H2O, and CO2 are found to increase with the ex- 
citation temperature (Lahuis & van Dishoeck 2000; Boonmaim 
et al. 2000). Detailed radiative transfer calculations will be 
needed to compare the results of our theoretical model with 
these observations. We note that abundance jumps similar to 
those predicted have been observed in the low-mass 'hot core' 
around IRAS 16293-2422 (Schoier et al. 2002). 

Our results show that a different chemistry operates in the 
surrounding cooler envelopes from that in the central hot re- 
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gions where icy mantles can be evaporated, as has been indi- 
cated (e.g., Millar et al. 1997; Doty et al. 2002). In order to 
see the quantitative differences comparable with observations, 
we calculate the molecular column densities off-centred from 
the clump centre by a distance ^off , making use of Eq. Q with 
the distance from the clump centre, r' - r - r^s replacing r. 
The results at f = 10"^ yr are listed in Table |2] Parameters of 
©beam = 0.25 pc and Toff = 0.4 pc are adopted for compar- 
ing the results with observations toward the halo associated 
with G34.3h-0.15 by Thompson et al. (1999), which are also 
included in Table |2 For comparison, the molecular column 
densities calculated around the clump centre by Eq. with 
t - 10"* yr are also listed together with the observations which 
are the same as those in Fig. 3. The results of the model calcu- 
lations are consistent with most of the observations. The ratio 
of the calculated molecular column densities around the clump 
centre to those off-centred range roughly from a factor 2, when 
the radial abundances of the species are almost constant, to one 
order of magnitude, in the case of parent or daughter species 
which are hardly generated by the gas-phase reactions without 
grain mantle evaporation. 

4.2.3. Detailed chemistry 

In this subsection, we discuss some of the hot core chem- 
istry, mainly focusing on the effects of the injection of parent 
molecules from grain mantles on the radial abundance profiles 
of daughter and other related species. 

Fig. 2(a) and Fig. 2(b) show the time evolution of column 
densities of the carbon- and oxygen-bearing species, while 
Fig. 4(a) and Fig.4(b) represent the radial profiles of their abun- 
dances at f - 10"* yr The results related to oxygen chemistry 
can be interpreted based on the detailed analysis by Charnley 
(1997). Water is the main parent molecule among oxygen- 
bearing species and mainly destroyed by HCO^ to produce 
H30^ in this model. This process keeps the HCO^ abundance 
low where H2O is injected from the grain mantles and is con- 
sistent with the result of Doty et al (2002). Hydroxyl molecules 
are formed by the dissociative recombination of HjO^, so that 
the OH abundance changes discontinuously at the H2O injec- 
tion radius, Hnj h,o- One of the destruction processes of car- 
bon monoxide is the reaction with He^, which makes atomic 
oxygen as well as atomic carbon. On the other hand, O and 
OH react with molecular hydrogen to reproduce water as O 

H2 Ht 

— > OH — ^ H2O at high temperature. Thus, the abundances 
of O and OH are suppressed in the inner hot region. These 
behaviours are also consistent with the results of Doty et al 
(2002). The reaction between OH and N generates NO, which 
makes the radial abundance profile of NO similar to that of OH. 
Molecular oxygen is not only injected from the grain mantles 
but also synthesized from O and OH. HCO is principally pro- 
duced by the destruction of H2CO by OH, so that the HCO 
abundance jump appears at rinj^H^co- C3H4 is chiefly formed 

HCO+ , CH4 , e 

mitiated from C2H2 as C2H2 C2H+ C3H+ C3H4 
in this model. Therefore, the C3H4 abundance increases over a 
radial distance between nnj^HjO^ inside which the HCO^ abun- 
dance is suppressed, and nnj^CjHz- 



Fig. 2(c) and Fig.4(c) are the time evolution of column den- 
sities and the radial abundance profiles of nitrogen-bearing 
species, whose detailed hot core chemistry have analysed by 
Millar et al. (1997) and Rodgers & Charnley (2001). The ni- 
trogen molecule is the most dominant nitrogen-bearing species 
in the initial gas in this model, and is converted into N and 
N^, through reaction with He^. Ammonia is one of the possi- 
ble nitrogen-bearing mantle species and is destroyed by ionized 
molecules to produce NH2. Meanwhile, it is formed from gas- 
phase reactions initiated from as — ^ NH^ — ^ NHj 

^ NH+ -— » NH+ ^ NH3. The radial abundance profile of 
HCN has two apparent discontinuous jumps in this model. One 
is at the radius where Tacc.HCN = Tevap.HCN and the other is 
at /"inj.HjCO as H2CO is one of the main precursors of CH2, 
which reacts with N to produce HCN. The major production 
reaction of HNC is that between NH2 and C. Thus, the HNC 
abundance jumps at nnj.NH,- CH3CN is mainly produced by the 
radiative association between HCN and CH3, so that the ra- 
dial abundance profile of CH3CN resembles that of HCN. The 
CH3CN abundance is suppressed inside nnj.HjO because H^, 
which helps produce CH^ from its reaction with CH3OH, de- 
creases in abundance as it preferentially transfers its proton to 
the abundant water molecule. HC3N is generated from N as N 

^ C3N ^ HC3N+ — » H2C3N+ HC3N. C3H is a product 
of C2H2, which causes the HC3N abundance jump at rinj,c,HT 
CN is also chiefly formed from H2C3N^ inside rinj.c,Hi, but the 
jump does not appear in the CN abundance because CN can be 
synthesized by neutral-neutral reactions such as Ch-NO outside 
'"inj.CHT The CN abundance is suppressed in the inner hot re- 
gion due to the loss by a reaction with H2 which creates HCN. 

Fig. 2(d) and Fig. 4(d) represent the time evolution of col- 
umn densities and the radial abundance profiles of sulfur- 
bearing species, respectively. Their basic behaviour can be 
explained following the analyses by Millar et al. (1997) and 
Charnley (1997). Hydrogen sulfide is a possible sulfur-bearing 
grain mantle species, and can be easily destroyed by atomic hy- 
drogen especially in the inner hot region because of the small 
activation energy barrier of this reaction. This drives the sulfur 

H H Oi OH 

chemistry in hot cores as H2S — > SH — > S — ^ SO — > SO2. 
Thus, the SO abundance is lowered inside rinj.H,o, where the 
OH abundance is high, while the SO2 abundance drops outside 
'"inj,H,o- Also, both the SO and SO2 abundances change dramat- 
ically at rinj,o,. NS is produced by the reaction between SH and 
N, so that the NS abundance changes discontinuously at h,s- 

c 

CS is mainly generated as SO — > CS in the inner region, while 
HCS+ CS in the outer region. On the other hand, HCS^ is 
synthesized by the proton transfer process to CS in the inner re- 
gion, whereas by the reaction between CS^ and H2 in the outer 
region. Therefore, the radial abundance profiles of HCS^ and 
CS resemble each other OCS is injected from grain mantles in 
this model in order to represent the observed high abundance of 
OCS molecular lines, and based on the observations of ice ab- 
sorption features by Palumbo et al. (1997) (see also discussions 
in Hatchell et al. 1998b). 

Fig. 2(e) and Fig. 4(e) show the time evolution of column 
densities and the radial abundance profiles of large oxygen- 
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Fig. 5. The same as Fig. 4 but with the trapping of mantle 
molecules in water ice. Large differences appear in the radial 
distributions of some molecules compared to those in Fig.4. 

Table 3. Initial abundances of mantle molecules for Model B 



Species Abundance 


Species 


Abundance 


Species 


Abundance 


C2H2 5.0 10-' 


CO2 


3.0 lO-** 


O2 


1.0 io-« 


CH4 2.0 10-' 


H2CO 


2.0 10-« 


H2O 


2.8 lO"'* 


C2H4 1.0 10-* 


CH3OH 


4.0 10-' 


NH3 


8.0 10-' 


C2H6 1.0 10-* 


C2H5OH 


1.0 IQ-* 


H2S 


3.0 IQ-' 








ocs 


1.0 IQ-^ 



r[pcl 



bearing species, which can be interpreted according to the 
analyses by Millar et al. (1991) and Charnley et al. (1992). 
Methanol is one of the possible large oxygen-bearing grain 
mantle molecules and is protonated by H^. The product of 
the protonation, CH3OHJ, leads to the formation of daugh- 
ter molecules as CH3OH2 
CH3OH 



H2C0 



HCOOCH+ HCOOCH3 



and CH^OHi 



C2H6OH+ 



(CH3)20. Thus, abun- 



dance jumps of HCOOCH3 and (CH3)20 appear at nnj^cH^OH- 
CH3CHO is mainly made by the reaction between atomic oxy- 
gen and C2H5 which is generated from C2H6, so that the 
CH3CHO abundance changes dramatically at nnj.CjH,,- CH2CO 



is chiefly formed from the reactions of C2H3 — > CH2CO and 
CO „^ e 



CH+ 



CH3CO"' —> CH2CO inside Hnj.c^Hi, while from the 
latter reaction otherwise. In this model C2H5OH is assumed to 
be created on dust grains in order to reproduce the observations 
of C2H5OH molecular Unes (cf. Millar et al. 1997, Rodgers & 
Charnley 2001, Ikeda et al. 2001). 

5. Discussion 

5.1. Trapping of mantle molecules in water ice 

It is known that water is one of the most abundant ice com- 
ponent in dense clouds. Some observations of ice absorption 
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Table 4. Molecular column densities of Model B and observa- 
tions. 



Species 






[cm-2] 




Species 




A'mol 


[cm -] 




Model 




Obs. 


Model B" 




Obs. 


CO 


2.0 


10^" 


> 


2.9 


10'^ 


H2S 


2.4 


10'^ 


> 


1.1 


10'" 


CjH 


3.2 


10'5 


> 


8.2 


10" 


SO 


9.7 


10'^ 


> 


2.4 


10'^ 


CH3C2H 


1.3 


10'2 




1.8 


10'^ 


S02 


2.3 


10'6 




9.3 


10'5 


H2O 


3.5 


10" 




1.3 


10'" 


cs 


1.5 


10" 




1.0 


10'" 


H2CO 


2.3 


10" 


> 


2.3 


10'^ 


ocs 


1.9 


10'" 


> 


1.8 


10'" 


HCO+ 


5.5 


10'^ 


> 


1.4 


10'5 


H2CS 


1.8 


10'* 


> 


3.8 


10'" 


HCO 


6.4 


10'3 


> 


5.3 


10'3 


HCS+ 


5.6 


10'2 


> 


3.0 


10'3 


NO 


1.6 


10'« 


> 


5.8 


10'5 


NS 


5.0 


10'^ 


> 


4.0 


10" 


NH3 


2.6 


lO'*' 


> 


2.7 


10'« 


CH3OH 


7.4 


10'" 




2.6 


10'" 


HNC 


1.3 


10'5 


< 


3.3 


10'* 


C2H5OH 


1.8 


10'5 




1.7 


10'5 


HCN 


5.0 


10'^ 




1.2 


10" 


CH2CO 


1.3 


10'* 




6.7 


10'* 


CN 


1.3 


10'^ 


> 


5.9 


10'" 


CH3CHO 


1.3 


10'* 




2.4 


10'* 


HC3N 


1.1 


10'^ 


> 


2.5 


10'3 


HCOOCH3 


1.0 


10" 




1.4 


10" 


CH3CN 


3.4 


10'2 


> 


2.4 


10'* 


CH3OCH3 


5.5 


10" 




3.7 


10'5 



" Calculated results at 10* yr. 
'' In the case of Obeam = 0.025 pc. 



features in infrared spectra towards young stellar objects have 
suggested that a certain amount of CO is frozen in H20-rich 
matrices (e.g., review by Allamandola et al. 1999), and labo- 
ratory experiments have shown that solid CO is able to diffuse 
into the porous structure of H2O and be trapped in it (Collings 
et al. 2003). These studies suggest that mantle molecules other 
than CO also could be trapped in water ice. Thus, in this sub- 
section we examine the effects of this trapping on the chemical 
structure of the hot core G34.3h-0.15. In order to see the ex- 
treme case, we put the injection radius of each mantle molecule 
as rinj,, = rinj.H,o, and assume that it is completely trapped in 
H2O matrices outside rinj,H,o, except CO, for which we assume 
that x(CO) - 0.3x(H2O) is trapped outside rinj,HjO based on 
observations (e.g., Chiar et al. 1998; cf. Doty et al. 2002). As 
a result of the time-dependent chemical calculation described 
in Sect. 3, we get the molecular column densities of the hot 
core consistent with most of the observations at around 10"* yr 
(Table 4). The initial condition listed in Table|5]is used in this 
calculation, which is slightly different from that in Table 1, in 
order to fit the results to the observations. Although there are 
few differences in the integrated molecular column densities, 
we find that very different structure appear in the radial pro- 
files of molecular abundances at around lO"* yr between the re- 
sults in Sect. 4 (Fig. 4: Model A) and those obtained from this 
calculation (Fig. 5: Model B). Especially the differences are re- 
markable in the following molecules: (1) the binding energies 
of the molecules or their parent species are low, that is, the in- 
jection radii listed in Table are large, such as H2CO, HCO, 
H2S, and SO, or (2) daughter molecules whose production are 
suppressed where H2O is evaporated from grain mantle, such 
as CH3C2H and CH3CN. For comparison, the resulting radial 
abundance profiles of SO and CH3CN from Model A (dotted 
lines) and B (solid lines) are plotted in the linear radial scales 
in Fig. 6. Detailed observations of radial profiles of molecular 




Model A 
Model B - 




Fig. 6. The comparison of the radial abundance profiles of (a) 
SO and (b) CH3CN in the linear radial scales between the mod- 
els without (dotted lines: Model A) and with (solid lines: Model 
B) the trapping of mantle molecules in water ice. We can find 
the difference in the location of jump in (a) and that the abun- 
dance is low through the clump in Model B in (b). 



abundances may make it possible to examine what fraction of 
molecules is trapped in water ice in cooler region of hot cores. 

5.2. Dependence of chemical structure on density 
profile - General model 

The comparison between the results of chemical calculation 
and the observations of various kinds of molecular lines will 
be useful to constrain the physical condition of massive star 
forming regions (e.g., Viti & Williams 1999). The density pro- 
file of hot molecular cores is one of the possible tracers of the 
evolution process just before and after massive star formation. 
In this section we will investigate the dependences of the chem- 
ical structure of hot cores on the density profile to discuss the 
physical condition of massive-star-forming clumps. 

The chemical structures of three different types of molecu- 
lar clupms are investigated here. Their density profiles are set 
by Eq. ([3 with core radii of rc - 0.01, 0.05, and 0.1 pc, and 
plotted in Fig. 7(a) with dotted, solid, and dashed lines, respec- 
tively. Making use of these density profiles, we calculate the 
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temperature profiles by solving the radiative transfer equation 
as described in Sect. 2.2. The resulting temperature distribu- 
tions are shown in Fig. 7(b). We can see from the figure that the 
temperature profiles directly reflect the difference in the density 
profiles. This is because the temperature profiles of hot cores 
are determined by the balance between the absorption of radi- 
ation from the central star and the reemission of the radiation 
by dust grains in the core as described by Eq. (|2ji. The inner 
regions of the cores are optically thick enough to the radiation 
from the central stars so that the photons from the stars diffuse 
via interactions with the surrounding dusty material. If we sim- 
ply assume that the dust opacity depends on the frequency of 
radiation as Ky oc v" and the density profile of a clump is set 
as p oc r^i\ for example, the temperature profile of the clump 
illuminated by a central star is derived as T oc 
the optically thick limit from the equation of local radiative 
equilibrium and the radiative transfer equation (e.g., Adams & 
Shu 1985). Meanwhile, since the evaporation process of grain 
mantle species is very sensitive to the temperature as shown in 
Sect. 3, it is to be expected that the density profiles affect the 
radial abundance profiles of molecules in hot cores via the tem- 
perature profiles. In order to see the dependence on the density 
models clearly, the temperature profiles in the region of grain 
mantle evaporation are plotted with linear scales in Fig. 7(c). 

The luminosity of the central star is another factor that in- 
fluences the temperature profile of the clump. As we can see 
from Fig. 7(b) and (c), the temperature distribution of the clump 



I \ \ 1 \ \ 1 

0.05 0.1 0.15 0.2 0.25 0.3 

r[pc] 

Fig. 7. (a) Density and (b) temperature profiles in the logarith- 
mic scales for the models with L, = lO^L© and rc - 0.01 (dot- 
ted lines), 0.05 (solid lines), and 0.1 pc (dashed lines). The 
temperature profiles directly reflect the difference in the den- 
sity profiles, (c) Temperature profiles in the linear scales repre- 
sent clearly the dependence on the density models around the 
temperatures of grain mantle evaporation. Profiles of the model 
with Lt - 10^ Lq and rc = 0.05 pc (dot-dashed lines) are also 
plotted for comparison. 

whose core radius is r^. - 0.05 pc and whose central star lumi- 
nosity is L, = IO^Lq (dot-dashed lines) is similar to that of the 
model with rc = 0.01 pc and = IO^Lq. The difference in 
luminosity will be, however, distinguishable by observing the 
spectral energy distribution. The hydrogen column density will 
not have very much effect on the temperature structure in the 
optically thick region because basically T(r) depends not on 
the absolute density but on the steepness of the density profile 
as we have shown in the above example. The optically thick 
region spreads (shrinks) if the hydrogen column density is high 
(low) so that it will affect the temperature profile at the outer 
radii. 

We calculate the gas-phase chemical reactions in the 
clumps with the above density and temperature profiles and 
L, = 10^ Lq by means of the chemical model in Sect. 3. 
Fig. 8 gives the resulting radial abundance profiles of some 
molecules at 10"^ yr in the linear radial scales (while those in 
Fig. 4 are scaled in the logarithmic manner). They are parent 
or daughter species and have dramatic changes at certain radii 
related to grain mantle evaporation. H2O, CO2, H2CO, H2S, 
C2H2, and CH3OH are possible parent molecules, and have 
abundance jumps at their injection radii. The daughter species 
HC3N, CH3CN, and HCOOCH3 are produced mainly from 
C2H2, HCN, and CH3OH, respectively, in this model, and their 
abundance jumps appear at the same radii as those of the par- 
ent species. Comparison between reasonably high resolution 
observations of abundance jumps and precise radiative transfer 
and physical-chemical models of hot cores may make it possi- 
ble to test the formation processes of these molecules. 

From Fig. 8 we can see that the abundance jump of each 
molecule occurs at a smaller radius in the clump with a smaller 
core radius rc, that is, with a steeper density and therefore a 
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Steeper temperature profile. The molecules in Fig. 8 can keep 
high fractional abundances only in the inner region of the 
clump (< O.lpc) if Tc - O.Olpc. This result suggests that the 
chemically rich hot cores, whose observational beam-averaged 
molecular column densities are large, have flat density profiles 
rather than steep ones. On the basis of the low and high mass 
star formation theories, the central density profiles of clumps 
are thought to be steep just after the star formation, at least if 
accretion of the surrounding material to the central protostel- 
lar cores is important (e.g., Shu 1977; Yorke & Bodenheimer 
1999; Masunaga & Inutsuka 2000; Yorke & Sonnhalter 2002). 
Subsequently, dissipation by outflows will let the clumps grav- 
itationally unbound and their central density profiles will be 
flattened (e.g., Shu et al. 1987; Nakano et al. 1995). In the case 
of massive star formation, powerful outflows with high mass 
loss rates as well as radiation pressure and strong winds from 
the stars may accelerate the dissipation process. If massive stars 
are formed as a result of coalescence of lower mass stars, the 
envelopes will have different density profiles from those of iso- 
lated star-forming clumps (e.g., Boimell et al. 1998). Statistical 
observations of various kinds of molecular lines towards hot 
cores could constrain their density profiles and the massive star 
formation theories. 

The effects of inflow motion of hot cores have not been 
considered in this paper. If the inflow has a velocity of a few 
km s~' as is estimated from the observed atomic and molec- 
ular lines (e.g., Keto 2002), it will not affect the structure on 
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Fig. 8. The radial abundance profiles of some parent and daugh- 
ter molecules of the models in the linear radial scales with 
L. = IO^Lq and = 0.01 (dotted lines), 0.05 (sohd lines), 
and 0.1 pc (dashed Unes). The abundance jumps appear at the 
inner radii in the models with the smaller core radii, that is, 
having steeper density and temperature profiles. 



the scale of about 0.25 pc, which corresponds to the observa- 
tional beam size for most of the molecular lines we have used in 
this paper, within the timescale of 10"* yr. Comparison between 
model calculations and the forthcoming high spatial resolution 
and high sensitivity observations could reveal the dependence 



H. Nomura and T. J. Millar: The physical and chemical structure of hot molecular cores 



13 



of the chemical structure of hot cores on the inflow motion. In 
addition, the temporal evolution of young massive stars in hot 
cores is also a factor that may possibly influence the fractional 
abundances of some species (cf. Viti & Williams 1999). 



6. Summary 

We have constructed the density and temperature profiles of 
hot cores that have embedded luminous objects, using radiative 
transfer calculations together with observations of the spectral 
energy distribution. We have used these to investigate hot core 
chemistry taking into account the temperature dependent grain 
mantle evaporation process. 

As a result of applying the model to the hot core 
G34.3+0.15 and calculating the time-dependent gas-phase re- 
actions in hot gas heated by a central massive star, it is con- 
firmed that the age of around 10^ yr is preferred to reproduce 
the large column densities of both the parent and the daughter 
molecules. In fact, the calculated column densities are consis- 
tent with most of those observed toward G34.3+0.15 around 
10"* years after central luminous star formation with a particu- 
lar initial composition for the grain mantle molecules. 

Also, we find that dramatic changes appear in the radial 
profiles of fractional abundances of some parent, daughter and 
related species at the parents' injection radii, inside which the 
evaporation time of the molecules from dust grains are shorter 
than their accretion time. The abundance jumps appear be- 
cause (1) most of the parents are difficult to destroy within 
the timescale of 10^ yr, and (2) the initial abundances of par- 
ents is large enough, so that despite their slow destruction 
rates inside the injection radii, they produce appreciable abun- 
dances of daughter products, while (3) these parent and daugh- 
ter molecules are hardly produced by the gas-phase reactions 
outside the injection radii. 

Finally, we discussed the influence of some chemical and 
physical properties on the chemical structure of hot cores. First, 
we investigated the efl'ects of the trapping of mantle molecules 
in water ice. As a result, we find large changes in the location of 
jumps in the radial abundance profiles of some species and that 
it changes completely the radial abundance profiles of daughter 
molecules whose production are suppressed in regions where 
H2O is evaporated from grain mantle. Also, we investigated the 
chemical structure of hot cores with different density profiles. 
Our results suggest that the chemically rich hot cores have flat 
density profiles rather than sharp ones. This is because a steeper 
radial density profile leads to a steeper temperature distribution 
as long as the clump is optically thick to radiation from the 
central star. This means that the abundance jumps appears at 
smaller radii and that therefore beam-averaged molecular col- 
umn density cannot be large. Observational possibilities of con- 
straining the density profiles of hot cores and the massive star 
formation theories are also suggested. 
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